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ABSTRACT 

Ultra-luminous X-Ray Objects are off-nucleus point sources with Lx = 10 39 - 
10 41 erg/s but the nature of such systems are largely unidentified. Here we report 
a 2.1 hour period observed in a Chandra ACIS observation for ULX M51 X-7, 
which is located on the edge of a young star cluster in the star forming region in 
a spiral arm. In two ACIS observations separated by one year, the ULX changed 
from a high-hard to a low-soft spectral state, in contrast to most Galactic low 
mass X-ray binaries. Based on its period and spectral behaviors, we suggest that 
this ULX is a low mass X-ray binary system, with a dwarf companion of 0.2-0.3 
M 0 and a compact accretor, either a neutron star or a black hole, whose mass is 
not well constrained. Relativistic beaming effects are likely involved to produce 
the observed high X-ray luminosities given its low accretion rate as inferred from 
a sustainable accretion scenario via Roche lobe overflow. 

Subject headings: Galaxy: individual(M51) — X-rays: binaries 


1. Introduction 

X-ray emission from active galactic nuclei and quasars, powered by accretion onto nu¬ 
clear super-massive black holes, are usually > 10 41 erg/s, while for the stellar binary systems 
in the Milky Way, the luminosities are typically 10 33 -10 38 erg/s. X-ray objects with lumi¬ 
nosities of 10 39 -10 41 erg/s have been observed in some external galaxies as off-nucleus point 
sources by X-ray satellites such as EINSTEIN, ROSAT, ASCA, and recently Chandra X-ray 
Observatory (e.g., Fabbiano 1989; Marston et al. 1995; Zezas et al. 2002). These objects 
have several names, a popular designation being “Ultra-luminous X-ray Objects”, or ULXs. 

One suggestion for the nature of these ULXs is that they are binary systems with 10 3 - 
10 4 M q black holes as primaries (Colbert and Mushotzky 1999). Such black holes, if they 
exist, are the missing links between stellar mass black holes and super-massive black holes 
in the nuclei of galaxies. This suggestion is consistent with the X-ray spectral of some 
ULXs (Makishima et al. 2000). However, the formation of such massive black holes is not 
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predicted by stellar evolution theory and it may be impossible to form such objects even in 
dense star clusters. Alternatively, these sources may be stellar mass black holes or neutron 
stars whose emission is beamed, thereby representing micro-quasars (King et al. 2001; King 
2002). If the emission is beamed, the intrinsic luminosities become sub-Eddington, and there 
are known examples of beamed Galactic X-ray sources (c.f., Mirabel & Rodriguez 1999). It 
is also possible that the luminosities are truly super-Eddington, for example, obtainable from 
accretion disks with radiation-driven inhomogeneities (Begelman, 2002). 

To identify which of the above suggestions is correct, one natural and fruitful way is to 
associate them with known classes of objects by comparing spectral behaviors, short-term 
and long-term variabilities, etc. Spectral analyses of ULXs, not without problems, generally 
support the mass-accreting black hole scenarios (e.g., Makishima et al. 2000). Spectral 
transitions were observed for two ULXs in IC 342 between high-soft and low-hard states, 
reminiscent of some black hole X-ray binaries in our Galaxy (Kubota et al. 2001). While 
many ULXs are variable during individual observations, few periods have been reported. 
So far, only a possible periodicity at several tens of hours of a LTLX in IC 342 has been 
reported (Sugiho et al. 2001). In this paper, we report a 2.1 hour period for a LILX in 
M51 (NGC 5194), M51 X-7 as designated in Roberts & Warwick 2000. This ULX is located 
on the outskirts of a young open cluster on a spiral arm of M51. In §2, we present the 
two Chandra ACIS observations utilized and the results. A period of 7620 seconds for the 
LILX is found from the X-ray light, and a spectral transition is found between two Chandra 
observations. In §3, we discuss the implications of its period and spectral evolution, and 
suggest that this LILX is a low mass X-ray binary system, with a dwarf companion of 0.2-0.3 
Mq and a compact accretor. Possible mechanisms to produce the high X-ray luminosity is 
also discussed. For the distance to M51, we use 7.7 Mpc (Roberts & Warwick 2000). 


2. Data analysis and results 

Two Chandra ACIS observations of M51, ObsID 354 (15 kilo-seconds, 2000-06-20) and 
ObsID 1622 (27 kilo-seconds, 2001-06-23), were retrieved from the archive. The data were 
analyzed with CIAO 2.2.1, CALDB 2.12, XSPEC 11.2.0 and XRONOS 5.19. Pixel random¬ 
ization of event positions was removed to recover the original positions before CIAO task 
wavdetect was used to detect discrete sources on ACIS-S chips. The accuracy of absolute 
positions for these sources is expected to be 076 (Aldcroft et al. 2000). 

The LILX M51 X-7 was observed in both data sets, with ACIS positions differing by 
(075, 071). For the true position, we take the average, i.e., RA=13:30:0.99, Dec=47:13:43.9, 
with a positional uncertainty of 074. Registering this position onto an HST WFPC2 image, 
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we find that it falls onto the edge of a star cluster in a dusty star forming knot on a spiral 
arm of M51. Within its positional error circle are a few stars fainter than 23 mag. The 
details of its potential optical counterpart will be reported in a future paper. 

For timing and spectral analyses, we selected circular apertures with such sizes that 
the number of source photons outside of the aperture is approximately the number of back¬ 
ground photons inside the aperture. The size chosen is 4"0 for Observation 354, and 2 / /0 
for Observation 1622. The light curve of X-7 in the first ACIS observation (ObsID 354), 
defined by about 550 photons, appears to have the structure of a single period, as shown 
in Figure 1 in contrast with the background light curve. When fitting a periodic signal, 
we derive a period of 7620 seconds (using XRONOS/efsearch with 10 second resolution). 
The uncertainty in the period is approximately 500 seconds (a more accurate value is not 
possible due to the limited number of periods involved). To determine whether this could be 
an instrumental effect, we carried out a number of tests. We formed the light curves of all 
other bright sources in the same ACIS held, but found no evidence for periods. The period 
cannot be caused by background variation, since the count rate of the background varies by 
less than 10% (and is not periodic), while the count rate of X-7 varies by an amplitude of 
more than 50%, for both soft photons (< 1.5 keV) and hard photons (> 1.5 keV). 

Dithering may cause spurious periods when a source is dithering across a CCD node 
boundary or gap, falls off of the detector, or crosses bad pixels periodically. In this observa¬ 
tion, the dither period is 1000 seconds in the Chip Y direction, and 700 seconds in the Chip 
X direction. The dither pattern is a Lissajous figure, spanning 16 arc-seconds from peak to 
peak. The ULX falls on node 2 of chip S3 at (556,459). It is 22 arc-seconds away from the 
node 1-2 boundary, and no bad pixels fall within 30 arc-seconds of the source. Thus it is 
unlikely that the period is due to the dithering period. 

In the second ACIS observation (ObsID 1622) that was taken 1 year after the first, only 
53 photons were collected for the ULX, too few to make a reliable periodicity analysis. Note 
that the second observation was longer than the first by a factor of about 1.8, so the photon 
flux in the 0.3-10keV decreased by a factor of about 18 between these two observations. 

Photons in the 0.3-8 keV band are used to perform spectral fits for both spectra, as shown 
in Figure 2. The CIAO thread acisspec was used to extract the spectra, and corrarf was 
used to correct the continuous degradation in the ACIS quantum efficiency. To get good y 2 
statistics, we grouped the first spectrum (ObsID 354) requiring a minimum of 20 photons per 
bin. For the second spectrum (ObsID 1622), we required a minimum of 5 photons per bin to 
get enough bins. We tried absorbed (wabs) power-law (powerlaw), absorbed multi-color disk 
(diskbb) and absorbed blaekbody (bb) models on each spectrum. In Table 1 we list the best 
fit parameters, the errors for the best-fit parameters are for the 90% confidence level. For 
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the first spectrum, an absorbed blackbody model is not acceptable, but an absorbed power- 
law model is acceptable, with a power-law photon index T = 1.251q 2 2 ; an( l an absorbed 
multi-color disk model is also acceptable, with an inner edge temperature kT in = 2.10to;4Q 
keV. To test whether there is an additional component due to a neutron star solid surface, 
we also fit an absorbed (multi-color disk + 2 keV blackbody) model to the first spectrum. 
With one more free parameter, this model improves the spectral fit marginally as compared 
to an absorbed multi-color disk model, and results in a smaller inner edge temperature, 
kT in = 0.82^q;3q keV. However, this additional component does not make a better model 
than a simple absorbed power-law model. 

For the second spectrum, an absorbed power-law model gives a power-law photon index 
of T > 2.6, an absorbed multi-color disk model gives an inner edge temperature of kT in = 
0.29 ± 0.09 keV, while an absorbed blackbody model is not acceptable. As compared with 
the first spectrum of higher photon flux, the second spectrum is much softer, and the best-fit 
absorbing column density is higher too. For simple absorbed power-law models, we switched 
the photon indices of the two spectra, and fitted a frozen V = 2.6 to the first spectrum, and 
a frozen T = 1.25 to the second spectrum, and found that both models can be accepted only 
on the levels < 10" 2 . Thus we conclude that the first spectrum is truly harder than the 
second spectrum. Absorbed power-law models with smaller frozen nn values were also fitted 
to the second spectrum and gave worse fits, indicating that the high column density in the 
second spectrum is real. 


3. Discussion 

Periods are not uncommon in Galactic X-ray point sources, such as cataclysmic variables 
(CVs), low mass X-ray binaries (LMXBs) and high mass X-ray binaries (HMXBs). HMXBs 
have periods of more than a few days, because they are usually widely separated systems. 
A two hour period is common for CVs, where the accretor and the donor are white dwarfs. 
However, CV systems have typical temperatures under 10 6 K, with most of the radiation 
emitted in the ultraviolet region and with X-ray luminosities below 10 34 erg/s. LMXBs 
usually have periods of 10+ hours, but they can have periods as short as one hour, for 
example, 1624-490 has a period of 0.69 hours. The source studied here, with a period of two 
hours and a temperature of ~ 1 keV, is most similar to the low mass X-ray binaries. 

The LMXB picture of LILXs does not contradict its location in a star forming re¬ 
gion. While observations show that LILXs preferentially occur in active star forming regions 
(e.g.,Zezas et al. 2002), indicating their short lifetimes, they nevertheless could be a short 
stage in the long lifetime of LMXBs. Indeed, some LILXs have also been found outside of 
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star forming regions, for example, in gas poor elliptical galaxies such as NGC 1399 (Angelini 
et al. 2001), and in a globular cluster on the outskirt of the bulge in edge-on spiral galaxy 
NGC 4565 (Wu et al. 2002). 

If X-7 is a LMXB with an orbital period of about 2 hours, there are a number of 
implications. The secondary star would be filling its Roche Lobe and undergoing mass loss 
to the primary, and from this picture, we can estimate the mass of the secondary, M\. 
For a mass ratio of q — Mi/M 2 < 0.8, the Roche Lobe radius is R cr = 0.46a( M A + M ) 1//3 
(Paczynski 1967), in which a is the separation between the donor and the accretor, and 
M 2 is the mass of the accretor. Combined with Kepler’s 3rd law, this leads to P or b = 
8.9(i?i/i? 0 ) 3 / 2 (M Q /M 1 ) 1 / 2 hours. For a late-type low mass star, if we adopt an analytical 
mass-radius relation Ri/R e = /M 0 (Verbunt 1993), we obtain the mass of the donor to 

be Mi ~ 0.23M 0 . We also tried to obtain the mass of the donor by adopting computed 
mass-radius relation from the Geneva stellar models (Lejeune et al. 2001). We extended the 
star mass down to 0.1 M 0 by extrapolating computed lgR-lgM relation for low mass stars 
(0.8 < M/Mq < 0.9). These computed stellar models give an equivalent analytical relation 
of Ri/Rq ~ ( Mi/Mq ) L2 , and the resulted donor mass is about 0.3 M 0 . Such a low mass 
secondary in M51 is impossible to identify even with the Hubble Space Telescope. 

The mass of the accretor cannot be constrained by the period under our assumptions. 
The accretor is most likely a neutron star, or a black hole , as a compact object is needed to 
produce the observed hard spectra. One can refine this suggestion further by comparing the 
observed spectrum with models of compact objects, as we try here. Neutron star LMXBs 
and black hole LMXBs at high luminosities might be distinguished by their spectral features 
(cf. Tanaka & Lewin 1995), since high luminosity neutron star LMXBs exhibit a multi-color 
disk component and a blackbody component with a color temperature of ~ 2 keV, which is 
interpreted as emission from the neutron star surface, while black hole LMXBs do not show 
such a single-color blackbody component at any luminosity level. As the luminosity decreases 
further, the spectrum becomes a single power-law, and distinction between black hole and 
neutron star LMXBs diminishes. For X-7, the model with an additional 2 keV thermal 
component improves the spectral fit marginally as compared to an absorbed multi-color disk 
model, which may be very weak evidence that the accretor is a neutron star. However, such 
a marginal improvement may be due to the additional freedom of parameters, and does not 
really reflect the preference of a thermal component. Also, its trend of spectral transition 
from a hard high state to a soft low one has some Galactic parallels, such as Nova Muscae 
1991 (Ebisawa et al. 1994), whose system is comprised of a low mass companion and a black 
hole of ~ 6 Mq. Such spectral transitions, as noted by Terishima & Wilson (2002), are in 
contradiction to most Galactic LMXBs which switch between soft high and hard low states. 
Thus it is premature to identify it as a neutron star or a black hole with current observations. 
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One concern about this picture is that mass transfer from the less massive donor to 
the accretor will widen the orbit and enlarge the Roche lobe with respect to the donor, 
and thus prevent further mass transfer via Roche lobe overflow. This can be avoided if 
the orbit angular momentum can be dissipated effectively by gravitational radiation, with 
the rate —( j / J)gr = + M 2 )/a 4 . The separation a evolves according to the 

equation a/a = 2 J/ J — 2(1 — M\ / M 2 ) M { / M l under conservative mass transfer via Roche 
lobe overflow(Verbunt 1993). In order for the orbit not to widen, the mass transfer rate via 
Roche lobe overflow must be less than a value of about ICR 9 — ICR 8 M@/yr for a primary of 
3-100 Mq. In a steady-state, accretion disks with these accretion rates generate luminosities 
of 10 3 ' —10 38 erg/s regardless of the mass of the accretor, less than 10 39 erg/s we observed for 
X-7. The larger observed energy would require an additional phenomenon, such as a transient 
outburst or relativistic beaming. For example, X-ray novae can have substantial outbursts, 
but usually with a duty cycle of less than 10%. However, this source has been brighter than 
10 39 erg/s during the past two decades (Figure 3), in conflict with expectations for dwarf 
novae. Therefore, the more likely explanation for the high luminosities is by relativistic 
beaming, as suggested by King et al. (2001). 

An alternative model for the binary and the periodicity is that the period arises due to 
a hot spot in the accretion disk around the primary. This picture removes the restrictions 
on the mass of the secondary and its separation from the primary. If we assume Keplerian 
motion of the hot spot around the accretor, this 2 hour period corresponds to a distance of 
~ 1 — 2 R 0 for an accretor mass of ~ 3 — 10 M 0 . At such distances, the multi-color disk 
temperature is normally 10 4 K. A hot spot at such distances fails to explain the fluctuation 
amplitude for hard photons (>1.5 keV). 

Another possibility is that the period is not real and results from the limited sampling of 
a randomly variable object. If the period reflects that of orbital motion, it should be highly 
reproducible, provided that the source is sufficiently bright. If it is due to the motion of a 
hot spot on a disk, one might expect the hot spot to fade with time (then the period would 
vanish while the source might still be luminous) or the spot might move inward in the disk, 
leading to a shorter period. Additional observations that could distinguish between these 
expectations would be greatly beneficial in clarifying the nature of the variation seen in M51 
X-7. 
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gratefully acknowledge support for this work from NASA under grants HST-GO-09073. 



- 7- 

REFERENCES 

Aldcroft, T.L., Karovska, M., et al. 2000, Proc. SPIE, 4012, 650 

Angelini, L., Loewenstein, M., and Mushotzky, R.F. 2001, ApJ, 557, L35 

Begelman, M.C., 2002, ApJ, 568, L97 

Colbert, E. J. M. and Mushotzky, R. F. 1999, ApJ, 519, 89 

Ebisawa, K., Ogawa, M. et al. 1994, PASJ, 46, 375 

Fabbiano, G. 1989, ARA&A, 27, 87 

King, A. R., Davies, M. B., Ward, M. J., Fabbiano, G. and Elvis, M. 2001, ApJ, 552, L109 
King, A. R., 2002, MNRAS, 335, L13 

Kubota, A., Mizuno, T., Makishima, K. et al. 2001, ApJ, 547, L119 

Lejeune, T., and Schaerer, D., 2001, A&A 366, 538 

Makishima, K., Kubota, A. Muzuno, T. et al. 2000, ApJ, 535, 632 

Mirabel, I. F. and Rodriguez, L.F., 1999, ARA&A, 37, 409 

Marston, A.P., Elmegreen, D., and Elmegreen, B. 1995, ApJ, 438,663 

Paczynski, B. 1967, Acta Astron. 17, 287 

Roberts, T. P. and Warwick, R. S. 2000, MNRAS, 315, 98 

Sugiho, M., Kotoku, J., Makishima, K., Kubota, A., Mizuno, T., Fukazawa, Y., and Tashiro, 
M. 2001, ApJ, 561, L73 

Tanaka, Y. and Lewin, W.H.G. 1995, in X-ray Binaries, eds. Lewin, W.H.G, van Paradijs, 
J. and van den Heuvel, E.P.J. P.265 

Terishima, Y. and Wilson, A. 2002, astro-ph/0204321 

Verbunt, F. 1993, ARA&A, 31, 93 

Wu, H., Xue, S.J., Xia, X.Y., Deng, Z.G., and Mao, S.D. 2002, astro-ph/0204359 
Zezas, A., Fabbiano, G., Rots, A.H., and Murray, S.S. 2002, astro-ph/0203175 


This preprint was prepared with the AAS LMRX macros v5.0. 



Table 1. Spectral Fits for M51 X-7 in 2 ACIS observations 


ObsID 

model 

N h 

(10 22 cm -2 ) 

r 

kT in 

(keV) 

kT 

(keV) 

xl/Aof 

absorbed Lx a 
(10 38 erg/s) 

unabsorbed L^ a 
(10 38 erg/s) 

354 

wabs* power law 

wabs*diskbb 

wabs*bb 

wabs* (bb+diskbb) 

0 06 +0 ' 07 
u.uo_o 06 

< 0.03 

0 02~*”°’ 09 
u - u ^_0.02 

1 25+ 0 - 24 
lzo -0.22 

9 19 +0.71 
ZliZ -0.42 

0-8 2 ti'.3O 

2.0 (frozen) 

0.95/21 

1.06/21 

3.47/21 

1.01/20 

18.6 

17.1 

18.8 

19.8 

17.1 

19.2 


wabs*powerlaw 

0 70" 1 " A ' U ^ 

u - ' u -0.55 

> 2.6 



0.82/5 

0.39 

9.34 

1622 

wabs*diskbb 

0 27"*" 1,09 


0.29 ± 0.09 


0.89/5 

0.30 

0.77 


wabs*bb 





4.61/7 




The X-ray luminosities are in 0.3-8 keV. 
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Light Curves: NGC 5194 X—7 & Background 

Bin time: 250.0 s 



Fig. 1.— Light curve of M51 X-7 in ACIS observation ObsID354. The ULX flux varies 
periodically by an amplitude of more than 50%. The overlaid light curve labeled by squares 
is the scaled background, which shows variation of less than 10%. Photons within 0.3-10 
keV were used for the light curves. 




































































Upper: ObsID 354 Lower: ObsJD 1622 



Fig. 2.— Spectral fits of M51 X-7 in two ACIS observations. In the first observation (Obsl 
354) the ULX was in a high hard state, while in the second observation (ObsID 1622, mod 
and y 2 in thick lines) it changed to a low soft state. 


channel energy (keV) 
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Days after first observation 

Fig. 3.— Long term variability of M51 X-7. Extracted from the archive were 1 EINSTEIN 
HRI observation (triangles), 1 ROSAT PSPC observation (pentagons), 5 ROSAT HRI obser¬ 
vations (squares) and 2 Chandra ACIS observations (circles). The instrumental count rates 
were converted to flux assuming absorbed power-law models (a) if in the high-hard state 
(open symbols), and (b) if in the low-soft state (filled symbols). A distance of 7.7 Mpc was 
assumed to convert fluxes to luminosities. 













